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!> ■ ABSTRACT 

We present first results on PLCKG266.6-27.3, a galaxy cluster candidate detected at a signal-to-noise ratio of 5 in the Planck All Sky survey. An 
CO , XMM-Newton validation observation has allowed us to confirm that the candidate is a bona fide galaxy cluster. With these X-ray data we measure 
an accurate redshift, z = 0.94 ± 0.02, and estimate the cluster mass to be M soa = (7.8 ± 0.8) x 10 14 M G . PLCKG266.6-27.3 is an exceptional 
\Q ' system: its luminosity of L x [0.5-2.0keV] = (1.4 ±0.05) x 10 45 erg s _1 equals that of the two most luminous known clusters in the z > 0.5 universe, 
and it is one of the most massive clusters at z ~ 1 . Moreover, unlike the majority of high-redshift clusters, PLCK G266.6-27.3 appears to be highly 
relaxed. This observation confirms Planck's capability of detecting high-redshift, high-mass clusters, and opens the way to the systematic study of 
population evolution in the exponential tail of the mass function. 

Key words. Cosmology: observations - Galaxies: cluster: general - Galaxies: clusters: intracluster medium - Cosmic background radiation, 
X-rays: galaxies: clusters 
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1 . Introduction decades, put strong effort into the discovery and characterisation 

of these objects. 

Very massive clusters above redshift z ~ 1, when the Universe 

was at half the present age, are predicted to be very rare. They Until recently it was possible to identify clusters of galax- 

potentially provide a sensitive probe to constrain deviations from ies only via optical/infrared or X-ray surveys. Indeed, the most 

the standard ACDM paradigm (e.g. Mortonson et al. 2011); e.g., distant clusters presently known have all been detected with 

owing to non-Gaussian perturbations, non-standard quintessence these techniques, e.g., the IR-selected cluster CLJ1449+0856 

models or modified gravity models (see Allen et al. 201 1, for a at z - 2.07 (Gobat et al. 201 1) and the X-ray selected system 

review). They are also ideal targets for studying key aspects of XMMU J105324.7+572348 at z=1.75 (Henry et al. 2010). For 

the gravitational physics that drives cluster formation, includ- both of these objects, extended X-ray emission has been detected 

ing measurement of the evolution of the mass concentration. For with XMM-Newton, confirming their status as fully established 

these reasons, the scientific community has, over the past two galaxy clusters; however, their total masses are more typical of 

systems in the poor cluster or group regime (< 10 14 M o ). Until 

* Corresponding author: M. Arnaud, monique.arnaud@cea.fr recently, the most massive cluster known in the z > 1 universe 
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Fig. 1 : Left panel: Planck isz map of PLCK G266.6-27.3 obtained with the modified internal linear combination algorithm (MILCA; Hurier et al. 
2010) with a spatial resolution of 10'. Middle panel: XMM-Newton exposure-corrected count rate image of the region indicated by the black box in 
the left panel. It is obtained using data from the EMOS 1&2 and EPN camera in the [0.3-2.0] keV energy band. The contours of the XMM-Newton 
image after wavelet filtering are overlaid in white. Right panel: corresponding XMM-Newton surface brightness profile. The green line indicates 
the best-fitting /?-model with a cusp (see text); the red line is this model convolved with the point spread function (PSF) of XMM-Newton, and the 
dashed line is the on-axis PSF of XMM-Newton, normalised to the central intensity. The source is clearly significantly extended. 



was XMMU J2235.3-2557 at z=1.39, discovered in the XMM- 
Newton Distant Cluster Project (XDCP) based on serendipi- 
tous cluster searches in XMM-Newton observations (Mullis et al. 
2005). For this system, Jee et al. (2009) estimate a mass of 
A^200 = (7.3 ± 1 .3) x 10 14 M from a weak lensing analysis. 

However, clusters are also detectable through the Sunyaev- 
Zeldovich (SZ) effect (Sunyaev & Zeldovich 1972), the spectral 
distortion of the cosmic microwave background (CMB) gener- 
ated via inverse Compton scattering of CMB photons by the 
hot electrons in the intra-cluster medium. Crucially, the total 
SZ signal is expected to be closely related to the cluster mass 
(e.g. da Silva et al. 2004), and its brightness insensitive to red- 
shift dimming. As a result, SZ surveys can potentially provide 
unbiased cluster samples that are as close as possible to be- 
ing mass-selected 1 . They offer an ideal way to identify massive, 
high-redshift clusters. One recent illustration is the detection of 
SPT-CLJ2106-5844 at z = 1.13 by the South Pole Telescope 
(SPT) survey (Foley et al. 2011). With an estimated mass of 
M 20 o = (1-27 + 0.21) x 10 15 M , SPT-CL J2106-5844 is nearly 
twice as massive as XMMU J2235.3-2557. 

The Planck 2 satellite has been surveying the sky in the mi- 
crowave band since August 2009 (Planck Collaboration 2011a) 
with a good (band-dependent) spatial resolution of 5 arcmin 
(Mennella et al. 201 1; Planck HFI Core Team 201 1). Compared 
to other SZ experiments such as ACT (Marriage et al. 2011) 
or SPT (Carlstrom et al. 2009), Planck brings unique nine-band 
coverage from 30 to 857 GHz and, most crucially, an exception- 
ally large survey volume. Planck is the first all-sky survey capa- 
ble of blindly detecting clusters (i.e., not guided in the search by 
prior observations), since the ROSAT All-Sky Survey (RASS) in 



1 In practice, the mass threshold detectable by Planck increases with 
redshift. The total SZ signal is not resolved by Planck at high z and it 
decreases with z due to the decreasing angular size of the object. 

2 Planck (http://www.esa.int/Planck) is a project of the European 
Space Agency (ESA) with instruments provided by two scientific con- 
sortia funded by ESA member states (in particular the lead countries: 
France and Italy) with contributions from NASA (USA), and telescope 
reflectors provided in a collaboration between ESA and a scientific con- 
sortium led and funded by Denmark. 



the X-ray domain. This coverage allows detection of the rarest 
clusters, the most massive objects lying in the exponential tail of 
the mass functitheron. 

Planck Collaboration (2011b) recently published the Early 
SZ (ESZ) sample, the first sample of galaxy clusters detected 
blindly in the all-sky maps from the first ten months of the 
Planck survey. The properties of this first sample already show 
that Planck is detecting previously unknown, massive clusters 
that do not appear in RASS or in other smaller area SZ surveys 
(Planck Collaboration 2011c). The ESZ comprises high signal- 
to-noise ratio (S/N > 6) Planck SZ sources up to z = 0.5. We re- 
port here on an SZ source that was blindly identified at S/N ~ 5 
in the Planck all-sky survey, and recent XMM-Newton validation 
observations confirm it is a massive cluster at z ~ 1 . 

In this paper, we adopt a ACDM cosmology with Ho = 
70 kms 1 Mpc" 1 , O m = 0.3, and Q A = 0.7. The factor E(z) = 
YOm(1 + z) 3 + £2a is the ratio of the Hubble constant at redshift 
z to its present-day value. The quantities Ms and Rg are the to- 
tal mass and radius corresponding to a total density contrast 6, as 
compared to p c (z), the critical density of the Universe at the clus- 
ter redshift; Ms = (4n/3) S p c (z) R 3 6 - The SZ flux is characterised 
by F500, where F500 D 2 A is the spherically integrated Compton pa- 
rameter within /?5oo (corresponding to S — 500), and Da is the 
angular-diameter distance to the cluster. 

2. Planck detection 

The blind search for clusters in Planck data relies on a multi- 
matched filter (MMF) approach (Melin et al. 2006). Candidates 
then undergo a validation process, including internal quality 
checks and cross-correlation with ancillary data and catalogues, 
as described in Planck Collaboration (2011b). This process 
produces a list of new Planck SZ cluster candidates above a 
given S/N threshold that require follow-up for confirmation. The 
XMM-Newton follow-up for validation, undertaken in Director's 
Discretionary Time via an agreement between the XMM-Newton 
and Planck Project Scientists, plays a central role in this confir- 
mation procedure. It consists of snapshot exposures (~ 10 ks), 
sufficient for unambiguous discrimination between clusters and 
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Fig. 2: XMM-Newton EMOS1 (black) and EMOS2 (red) spectra ex- 
tracted from a circular region of 1(5 in radius and centred in the cluster 
X-ray peak. Only the data points above 2 keV are shown for clarity, but 
data down to 0.3 keV are used in the spectral fitting. The line is the ther- 
mal model for the best-fitting redshift, z = 0.94 ± 0.02. The position of 
the redshifted Fe K line is marked. 



false candidates (Planck Collaboration 201 lc). The results of the 
first two runs (completed in September 2010) are reported by 
Planck Collaboration (201 lb,c). 

The XMM-Newton validation programme is continu- 
ing to explore lower S/N and detection quality criteria. 
PLCKG266.6-27.3, detected at S/N = 5.03, was observed in 
the framework of the third run of the XMM-Newton valida- 
tion programme, for which the analysis is on-going. This run 
comprises a total of 11 candidates detected at 4.5 < S/N < 
5.3 from the same Planck HFI maps used for the construc- 
tion the ESZ sample. The 1 1 candidates were sent for schedul- 
ing in November 2010 and the observations were performed 
between 22 December 2010 and 16 May 2011. Interestingly, 
PLCK G266.6-27.3 has been independently detected in the SPT 
survey. Its Planck position (6 h 16 m 6 s :6, -57° 47' 29") is consis- 
tent with that of SPT-CLJ06 15-5746 (Williamson et al. 2011, 
published on arXiv.org in January 2011, with a photometric 
redshift of z p h t = 1 +0.1). 



3. XMM-Newton validation 

3.1. Observation and data reduction 

PLCKG266.6-27.3 was observed with the XMM-Newton EPIC 
instrument (Turner et al. 2001; Striider et al. 2001), using the 
thin filters and the extended full frame mode for the 'pn-CCD' 
camera. The data analysis and validation procedure is described 
in Planck Collaboration (201 lc). Calibrated event lists were pro- 
duced with vl 1 .0 of the XMM-Newton Science Analysis System. 
Data that are affected by periods of high background due to soft 
proton flares were omitted from the analysis, and the remain- 
ing data were PATTERN-selected and corrected for vignetting, as 
described in Pratt et al. (2007). Bright point sources were ex- 
cised from the data. Background treatment is described in Pratt 
et al. (2010). In the spectroscopic analysis, the cluster compo- 
nent was modelled with an absorbed thermal emission model 
(mekal) with a hydrogen column density fixed at the 2 1 -cm value 
of Dickey & Lockman (1990). 



The observation, OBSID = 0658200101, was affected by 
soft proton flares. The net exposure time after flare cleaning is 
only 2.4 ks for the pn-CCD camera, with a particle background 
30% higher than nominal. The MOS camera data are less af- 
fected with a clean time of ~ 12 ks and a background excess 
about two times lower. We undertook a conservative approach to 
analysing spectroscopic data, since they are the most sensitive to 
the background estimate. We first fitted the data from the three 
cameras simultaneously, then fitted only the MOS cameras. The 
uncertainties in the physical quantities below reflect the differ- 
ence in best-fitting values between the two analyses and their 
errors. 



3.2. Confirmation and z estimate 

In Fig. 1 we show the vignetting-corrected count rate image of 
the cluster in the [0.3-2.0] keV band. An extended X-ray source 
is clearly coincident with PLCKG266.6-27.3. Its total EPIC 
count rate in the [0.3-2.0] keV band is (0.52 + 0.02)count/s 
within 2.3', the maximum radius of detection. The offset be- 
tween the X-ray cluster centre, defined as the emission peak at 
6 h 15 m 51 s :7, -57° 46'52'.'8, and the Planckblind position is 2'.07, 
consistent with the position reconstruction uncertainty, driven 
by the Planck spatial resolution and the source S/N (Planck 
Collaboration 201 lb). The extended nature of the source is con- 
firmed by comparing the surface brightness profile with the 
XMM-Newton point spread function (PSF) (Fig. 1, right panel). 
A typical (PSF-convolved) cluster surface brightness model con- 
sisting of a j6-model with a central cusp (Eq. 2 in Pratt & Arnaud 
2002) provides a good fit to the data and further supports the ex- 
tended nature of the source (Fig. 1). 

We extracted a spectrum within a circular region correspond- 
ing to the maximum significance of the X-ray detection (6 < 
1.'5). The iron K line complex is clearly detected (Fig. 2). Its sig- 
nificance is 3.6cr, estimated from a fit of the spectrum in the [2- 
6] keV band with a continuum plus a Gaussian line model. Since 
the centroid of the line complex depends on the temperature, 
the redshift is determined from a thermal model fit to the full 
spectrum, as described in detail in Planck Collaboration (201 lc). 
This yields a precise redshift estimate of z = 0.94 + 0.02. 



4. Physical cluster properties 

4.1. An exceptionally luminous and massive cluster 

We derived the deprojected, PSF-corrected gas density profile 
from the surface brightness profile, using the non-parametric 
method described in Croston et al. (2006). Global clus- 
ter parameters were then estimated self-consistently within 
^500 via iteration about the Msoo-^x relation of Arnaud 
et al. (2010), assuming standard evolution, E(z) 2 ^Msm = 
10 i4.567 ± o.oio [_^_j a561±0 018 Mq The quantity Fxj intro . 

duced by Kravtsov et al. (2006), is defined as the product of 
M g ^oo, the gas mass within /?5oo, and Tx- 7x is the spectroscopic 
temperature measured in the [0. 15-0.75] R500 aperture. In addi- 
tion, L500, the X-ray luminosity inside R500, was calculated as 
described in Pratt et al. (2009). All resulting X-ray properties, 
including iron abundance, are summarised in Table 1 . 

PLCK G266.6-27.3 is an exceptionally luminous system. Its 
[0.1-2.4] keV band luminosity of (22.7 + 0.8) x lO^ergs" 1 is 
equal to that of the fifth most luminous object in the MCXC 
compilation of Piffaretti et al. (2010), MACS J07 17.5+3745 at 
z = 0.55, discovered in the RASS by Edge et al. (2003). 
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Table 1 : Physical properties of PLCK G266.6-27.3 derived from 
XMM-Newton data. 



Parameter Value 



z 


0.94 + 0.02 


Abundance 


U.44 + U. 1 / solar 


^500 


0.98 + 0.03 Mpc 


M 500 


7.8+™ x 10 14 M o 


Yx 


1.10 + ° y° x 10 15 MokeV 


T x 


10.5 + [*keV 


t(<r 50 o) 


11.4 + |' 4 keV 


L 50 o([0.5-2.0]keV) 


14.2 + 0.5 x 10 44 ergs" 1 


L 50 o([0. 1-2.4] keV) 


22.7 + 0.8 x 10 44 ergs-' 



Notes. The M 500 -Y x relation with self-similar evolution is used to esti- 
mate R 5m (see Sect. 4). 

Moreover, its [0.5-2.0] keV band luminosity is consistent with 
that of SPT-CL J2106-5844, the most luminous cluster known 
beyond z — 1 (Foley et al. 2011). Collectively, these three clus- 
ters are the most luminous systems at z > 0.5. They are only 40% 
fainter than RXJ 1347.5-1145, the most X-ray luminous cluster 
known in the Universe (Piffaretti et al. 2010). 

Consistent with expectations for high-redshift Planck- 
detected clusters, we find that this cluster is extremely hot, Tx ~ 
llkeV, and massive, with a mass of M 50 o = 7.8^ x 10 14 M o . 
Our mass estimate is consistent with the less precise value, 
M 5 oo = 8 + 2 (statistical) + 1 .9 (systematic) x 10 14 M o , which is 
derived by Williamson et al. (2011) using the relation between 
SPT S/N and mass. Comparison of the masses of high-redshift 
systems is not trivial, because the estimation strongly depends 
on method, e.g which mass proxy is used and at what reference 
radius the mass is measured. On the basis of the published mass 
estimates, PLCK G266. 6-27. 3 would appear to be the most mas- 
sive cluster at z ~ 1 . Using the same factor to convert M500 to 
M 2 oo as Foley et al. (2011), we obtain M 200 = 15.5^xl0 14 M G , 
to be compared to M 200 = (12.7 ± 2.1) 3 x 10 14 M o for SPT- 
CL J2106-5844. However, the last value was derived by combin- 
ing X-ray and SZ data. A more direct comparison of M500 values 
estimated from the Msqq-Yx relation indicates that they are iden- 
tical within their uncertainties: M 5m = (7.8 + 0.8) x 10 14 M Q for 
PLCKG266.6-27.3 and M 500 = (9.3 + 2.0) x 1O 14 M for SPT- 
CL J2106-5844. 

4.2. F500 Compton parameter versus Y x 

The MMF blind detection was performed using the universal 
pressure profile of Arnaud et al. (2010) as a spatial template, 
leaving the position, characteristic size, 6500, and SZ flux, F500, 
as free parameters. The resulting flux, F500 = (5.6 + 3.0) x 
10~ 4 arcmin 2 , is consistent with the value, F^L = 6.4^[ q x 

10~ 4 arcmin 2 , expected from the measured value of Fx using 
the scaling relation derived from the universal pressure pro- 
file (Arnaud et al. 2010, Eq. 19). The cluster size, compara- 
ble to Planck's spatial resolution, is poorly constrained, #500 = 
3'.3 ±2^8. As discussed in Planck Collaboration (201 lb), the un- 
certainty on the blind F500 value is then large because of the 
flux-size degeneracy, where an overestimate of the cluster size 
induces an overestimate of the SZ signal. The SZ photometry 

3 The error includes an extra ~ 15% error accounting for uncertain- 
ties in the scaling relations. 



Table 2: SZ flux derived from Planck data with the reference value 
indicated in boldface. 



Method 


Definition 


Value 


#500 






(10~ 4 arcmin 2 ) 


(arcmin) 


MMF blind 


*500 


5.6 ±3.0 


3.3 ±2.8 


PWS blind 


F500 


6.5 ± 1.8 


3.9 ± 1.6 


MMF X-ray prior 


^500 


4.1 ±0.9 


fixed 


PWS X-ray prior 


F50O 


5.3 ±0.9 


fixed 


MILCA 


*tot 


5.9 ± 1.0 





Notes. Uncertainties on the blind values take into account the size un- 
certainty. 



can be improved by using the more precise XMM-Newton po- 
sition and size in the flux extraction. The F500 value obtained 
using these X-ray priors, F500 = (4.1 ± 0.9) x 10 -4 arcmin 2 , is 
lower than the value expected from the X-ray data at the 1 .7<x 
significance level. 

To check the robustness of the F500 estimate, we compared 
the MMF value with the one derived from the PowellSnakes 
(PWS; Carvalho et al. 2009, 2011) algorithm and the modi- 
fied internal linear combination algorithm (MILCA; Hurier et al. 
2010). The values are given in Table 2. PWS is a blind detec- 
tion algorithm that assumes the same profile shape as MMF, but 
is based on a Bayesian statistical approach, as fully described 
in Carvalho et al. (2009). MMF and PWS give consistent re- 
sults, the difference between MMF and PWS F500 values be- 
ing about 1.3 times the respective lcr uncertainties. MILCA is a 
component separation method that allows reconstruction of the 
SZ map around the cluster from an optimised linear combination 
of Planck HFI maps. In contrast to the MMF and PWS methods, 
the SZ flux derived from MILCA is obtained from aperture pho- 
tometry, i.e., with no assumptions on SZ profile shape or size. 
Assuming a typical conversion factor of 2/3 based on the univer- 
sal profile to convert the total F tot MILCA measurement to F500, 
the MMF and MILCA estimates are in excellent agreement. 

Several factors may affect the X-ray and SZ flux measure- 
ments and bring them out of accord. We have checked for pos- 
sible AGN contamination that could lower the Fgoo value us- 
ing the NVSS (at 1.4 GHz, Condon et al. 1998) and SUMSS 
(at 0.84 GHz, Bock et al. 1999) catalogues, but no bright radio 
sources are found in the cluster vicinity. The closest radio source 
with significant flux density is at 12'6 away. The source has a 
1 GHz flux density of 0.46 Jy. We also find no evidence of ra- 
dio contamination in the low-frequency Planck bands. On the 
other hand, the Fx measurement may also be increased by AGN 
contamination, from cluster members or foreground/background 
galaxies. Point source contamination is difficult to estimate ow- 
ing to the XMM-Newton PSF So, we estimate a maximum con- 
tribution to the X-ray luminosity from a central active galaxy of 
< 20%, assuming a point source model normalised to the central 
value of the X-ray surface brightness. The contribution to the 
gas mass, hence to Fx, would be less, provided that the source 
is not hard enough to significantly affect Tx- Nevertheless, only 
high-resolution X-ray imaging (e.g., from Chandra) can defini- 
tively establish whether X-ray AGNs at the cluster location af- 
fect our luminosity or mass measurement. A departure from the 
universal pressure profile would change the F500/FX ratio. The 
density profile shown below does not show any indication of 
this effect; however, deep spatially resolved XMM-Newton and 
Chandra spectroscopic observations are needed to derive the ra- 
dial pressure gradient from the core to R500. A final interesting 
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Fig. 3: Scaled density profiles of the REXCESS local cluster sample 
(Bohringer et al. 2007; Croston et al. 2008). The blue lines show the 
profiles for the cool core systems, and the orange lines the density pro- 
files of the morphologically disturbed systems. The density profile of 
PLCKG266.6-27.3 is shown with a thick red line. 



possibility is that gas clumping could affect the Fx measure- 
ments. A combination of X-ray and higher resolution SZ ob- 
servations is required to assess this point. 

4.3. Dynamical state and self-similarity of shape up to high z 

The available information indicates that PLCKG266.6-27.3 
may be particularly dynamically relaxed. The cluster image 
(Fig. 1, middle panel) does not show any sign of disturbance: the 
surface brightness is quite regular and quasi-azimuthally sym- 
metric within Rsoo- The offset between the X-ray surface bright- 
ness peak and the cluster brightest galaxy (Williamson et al. 
2011,Fig.l9) is less than 5". 

To further examine the dynamical state of the cluster, in 
Fig. 3 we compare its scaled density profile to those of clus- 
ters in the local representative X-ray-selected sample REXCESS 
(Bohringer et al. 2007; Croston et al. 2008). The radii are scaled 
by R500 and the density by the mean within R500. As extensively 
discussed by Pratt et al. (2009) and Arnaud et al. (2010), mor- 
phologically disturbed (i.e., merging) systems have systemati- 
cally shallower density profiles than more relaxed cool core ob- 
jects. This is illustrated in Fig. 3, where we indicate the scaled 
density profiles of the more relaxed cool core and the dynam- 
ically active merging clusters. The scaled density profile of 
PLCKG266.6-27.3 lies between the two classes, but with an 
indication of being closer to the relaxed rather than the merg- 
ing systems. It is thus possible that PLCK G266.6-27.3 is a cool 
core object at z ~ 1. Such objects are expected to be rare (e.g., 
Vikhlinin et al. 2007; Santos et al. 2010), and no cluster at this 
redshift has yet been found to contain a resolved central tem- 
perature drop that would confirm the presence of a cool core. A 
deep exposure at Chandra spatial resolution is needed to check 
this hypothesis. 

It is worth emphasising the similarity beyond the core of the 
density profile of this cluster with respect to REXCESS systems. 



This is the first piece of evidence for a similarity of shape up to 
redshifts as high as z ~ 1 . 

5. Conclusion 

PLCK G266.6-27.3 is the first blindly discovered Planck cluster 
of galaxies at z ~ 1. It has been confirmed by XMM-Newton 
in the framework of the on-going validation DDT observations. 
XMM-Newton data allowed us to measure the redshift with high 
accuracy (z = 0.94 + 0.02) and estimate the cluster mass to be 
M 500 = (7.8 + 0.8) x 10 14 M o . This XMM-Newton confirmation 
and redshift estimate is a clear demonstration of the capability 
of Planck for detecting high-z, high-mass clusters. 

PLCKG266.6-27.3 is an exceptional system, both in terms 
of its luminosity and its estimated mass. Furthermore, unlike 
other high-redshift clusters, it is likely to be a relaxed system, 
potentially allowing accurate hydrostatic mass measurements. It 
is thus a perfect target for deep multi-wavelength follow-up to 
address such important cosmological issues as the evolution of 
dark matter profiles, the evolution of the mass-Fsz relation, gas 
clumping, and the bias between X-ray and lensing mass esti- 
mates at such high redshift. 
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